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A New  Adaptive  Observer 


Kumpatl  S.  Narendra  and  Lena  S.  Valavani 


l\  introduction: 


An  adaptive  obaerver  in  defined  as  one  which  estimates  the  state  variables 

and  parameters  of  an  unknown  stable  linear  time-invariant  plant  from  its  input- 

output  data^  During  the  period  1973-1975  several  seemingly  different  versions  of 

jthe  adaptive  obaerver  [1-3]  appeared  in  the  control  literature.  Recently, 

| however,  it  has  been  shown  [5]  that  all  these  results  can  be  derived  in  a unified 
' — 

manner.  At  the  preaent  tin* there  are  two  distinct  approaches  to  the  design  of 
adaptive  observers  for  a plant  whose  input  output  behavior  can  be  represented  by 

h-^K.°r^*r  ®9u®t^on*  In  the  first  approach .,the  observer  is  of  the 

aaae  order  as  the  plant  and  is  referred  to  as  a minimal  (order)  observer.  Using 
the  second  approach, a non-mlnlmal  obaerver  of  order  (2n-l)  is  obtained.  Minimal 
observers  are  considerably  more  difficult  to  synthesise  than  non-mlnlmal  observers 
and  require  the  generation  of  additional  signals  for  the  stabilization  of  the 
adaptive  loop^J^T.  However,  they  have  the  advantage  of  yielding  simultaneously 


both  parameter  and  state  es time tea  of  the  plant.  Non-mininal  observers  are  con- 
siderably simpler  in  structure  but  the  n state  variables  of  the  plant  have  to  be 

Mtlaitid  from  the  available  (2n-l)  state  variables  of  the  observer. 

,t  > *6  prop oseef 

In  this  brief  paper,  we  peep see \a  new  observer7' which  appears  to  combine  the 


advantages  of  the  two  types  of  observers  described  above.  With  this  observer, the 
parameter  estimates  of  the  plant  are  directly  obtained  with  a structure  which  is 
no  more  complex  than  that  of  a non-mlnlmal  observer ^yf'which  is  widely  used  at 


the  present  time.  The  parameter  estimates  are  simultaneously  used  to  determine 
dirtetly  the  state  estimates  of  the  plant.  Under  certain  condltions^dAsewesed— 4 
eesedmm-S;  the  new  observer  has  a faster  rate  of  convergence  than  the  observers 
known  at  present, which  makes  It  particularly  attractive  for  use  in  the  control 
problem  JgrfC 
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From  a theoretical  point  of  view, the  lemmas  In  section  4 contain  the  prin- 
cipal ldeae  In  this  paper.  A direct  application  of  the  lemmas  results  In  a new 
prototype  for  an  error  model  and  this  In  turn  forms  the  mathematical  basis  for 
the  adaptive  observer. 

2.  The  Problem: 

A plant  P has  an  admissible  Input  u(t)  (a  scalar  valued  piecewise  continuous 

T 

function}  and  a corresponding  output  y(t).  A linear  system  {h  ,A  ,b  } is  said 

P p 

to  characterize  this  plant  If  every  pair  (u(t),y(t)}  satisfies  the  equations 

x - A x (t)  + b u(t) 

P P P P 

t (1) 

y(t)  - h x (t) 
p 


where  A Is  an  (rum)  matrix  and  h and  b are  constant  n-vectors.  The  identifica- 
P P 

cion  problem  may  be  qualitatively  defined  as  the  problem  of  constructing  a suit- 
able model  which,  when  subjected  to  the  same  Input  u(t)  as  the  plant, produces 

an  output  y (t)  which  tends  towards  y (t)  asymptotically. 

a P 

In  the  simplest  case,we  are  Interested  in  output  Identification  In  which 
|«1(t)|  ♦ 0,  where  e^t)  - y(t)  - ya(t),  for  * givon  Input  u(t).  When  the  input 
Is  "sufficiently  rich"  this  also  Implies  transfer  function  Identification.  If 
the  transfer  function  of  the  plant  Identified  corresponds  to  a unique  parametri- 
ze t Ion  of  the  plant, we  also  have  parameter  Identification. 

We  are  Interested  hers  In  a specific  parametrlzatlon  of  the  plant  and,hence. 
In  both  the  output  and  parameter  Identification  problems. 

The  plant  Is  described  by  equation  (1)  where  A^  la  a stable  (nxn)  matrix 
which  Is  In  output  canonical  form  so  that: 

■■  • ['K-] 


I 


and  h - [1,0,0, ... ,0) . The  vectors  a and  b In  aquation  (1)  are  the  unknown 

P P 

paraaeters  of  the  plant  which  have  to  be  estimated. 

If  Ab  Is  a known  stable  aatrlx  In  output  canonical  form  such  that 


' [*■  |— 


equation  (1)  may  be  expressed  as 


x ■ A x + 0.x,  + b u 

p ■ p 1 lp  p 

T 

x,  • h x and  0,  - a - a 
lp  p 1 p B 


To  Identify  the  paraaeter  vectors  a and  b , we  set  up  an  observer  which  Is 

P P 

described  by  the  differential  equation 

x - Ax  + bu  + d[0?v*  + ©Tv2] 

a a a a 1 2 

T ^ 

x,_  - h x 


where  and  d are  known  constant  vac tors,  0^(t)  and  0^(t)  are  n-dlaenslonal  ad- 

Jua table  paraaeter  vectors  which  are  the  estlaates  of  the  paraaeter  error  vectors 

1 2 

0^  - and  ©2  - (b^-b^) , and  v (t)  and  v (t)  are  known  Input  signals  which 

can  be  generated  using  the  Input  and  output  of  the  plant. 

If  x (t)  - x (t)  - e(t),  the  state  error  vector  and  hTe(t)  - e.(t),  the 

pa  l 

output  error,  we  have  froa  equations  (4)  and  (5) 

* “ V + Vl  + ®2U  ‘ dI®lyl  + ®2y2) 

- (6) 

ex  - h e 


The  problaa  now  la  to  determine  how  0^(t)  and  02(t)  aust  be  adjusted  so  that 
they  evolve  to  0^  and  02  respectively  and  result  In  the  output  error  e^(t)  tending 
to  aero  as  t ♦ •.  As  In  all  adaptive  observers, the  adaptive  lawe  are  expressed 


In  tim  of  the  tins  derivatives  of  8^(t)  end  e2(t). 

3.  Error  Equations  for  R»«imel  «nd  Non-mlnlmal  Adaptive  Obaervers : 

Equation  (6)  describes  the  behaviour  of  the  output  error  between  plant  and 
modal  In  terns  of  the  unknown  parameter  error  vectors  6.  and  6.  and  the  adjustable 

^ ^ e e 

A A A A 

parameters,  9^(t)  and  e2(t).  the  adaptive  lows  for  updating  8^(0  and  82<t) 

ere  specified,  equation  (6)  can  be  written  entirely  In  terns  of  the  output  and 
parameter  errors  of  the  overall  systen.  The  analysis  of  the  stability  behaviour 
of  the  error  equations  Is  obviously  central  to  the  understanding  of  the  conver- 
gence properties  of  the  adaptive  observer.  Asymptotic  stability  of  the  error 
equations,  for  example,  implies  that  the  states  and  parameters  of  the  observer 
converge  to  those  of  the  plant  asymptotically. 

Before  proceeding  to  consider  the  error  equations  of  the  new  adaptive  ob- 
server, we  shall  discuss  briefly  In  this  section  the  error  equations  which  have 
arisen  In  the  past  In  the  design  of  minimal  and  non-mlnlmal  observers. 

As  pointed  out  by  Narendra  and  Kudva  [5]  and  Auderson  [7],  an  Interesting 
prototype  for  error  equations  has  the  form  [Figure  1] : 


c • A^e  + d8  (t)v(t) 
T 

- h e 


(7) 


Where  e(t)  la  an  n-vector  of  errors,  6(t)  Is  an  m-vector  of  parameter  errors,  v(t) 
la  an  m-vector  of  piecewise  continuous  bounded  time  functions  which  ere  linearly 
Independent  on  the  semi- Infinite  Interval,  h and  d are  constant  n-vectors,  A^  Is 
an  (nxn)  stable  matrix  end  the  triple  {hT,Aa,d}  Is  completely  controllable  and 
completely  observable. 

T -1 

Zt  has  been  shown  that  If  the  transfer  function  h (sI-A^)  d Is  strictly 
positive  real  than  the  adaptive  law 


8 ■ -TtjV 


r - r > o 


(8) 





Figure 


■■■■■[ 
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vhare  r is  any  symmetric  positive  definite  matrix, will  result  in  a stable  system 
in  which  c(t)  will  tend  to  zero  and  6(t)  to  some  constant  vector  0q  as  t -*■ 

The  adaptive  lav  (8)  is  obtained  by  selecting  a function 

V(e,8)  - y[eTPe  + (9) 

* * dV 

as  a candidate  for  a Lyapunov  function  and  choosing  6(t)  to  make  V(e,8)  - — 
negative  sealdeflnlte.  It  haa  further  been  shown  that,  if  the  Input  v(t)  is 
"sufficiently  rich", the  parameter  error  vector  0(t)  will  tend  to  zero  [8]. 

The  above  simple  but  powerful  result  can  be  directly  applied  to  the  design 
problem  of  non-mininal  observers.  In  the  discussion  that  follows  It  is  shown 
that  It  can  also  be  related  to  an  error  differential  equation  (10)  which  arises 
In  the  design  of  minimal  observers. 

In  the  differential  equation 

«(t)  - Ame(t)  + 0(t)z1(t)  ei(t)  - hTe(t)  (10) 

let  e(t)  be  an  n-vector,  A an  (nxn)  stable  constant  matrix,  0(t)  an  n-vector 

n 

of  parameter  errors  and  z^(t)  a uniformly  bounded  scalar  function  of  time.  It 
is  desired  to  adjust  8(t)  using  only  available  input  and  output  data  (l.e.  e^(t) 
and  z1(t))  such  that  e^(t)  will  tend  to  zero  asymptotically. 

The  form  of  equation  (10)  does  not  lend  Itself  easily  to  the  generation  of 
adaptive  laws  of  the  fora 

e - f(a1,z1).  (11) 

However,  it  has  besn  shown  [4]  that  vectors  w(t)  and  v(t)  exist  such  that 

i - Aaa  + 8(t)s1(t)  + v(t)  ex(t)  - hTe(t>  (12) 

i - kut  + deT(t)v(t)  ei(t)  - hTEl(t)  (13) 
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have  the  same  outputs  e^(t)  and  e^(t).  The  vector  v(t)  can  be  generated  by  a 
dynamical  system  whose  Input  is  z^(t)  while  w(t)  depends  on  both  e^(t)  and  v(t). 

In  view  of  the  Input-output  equivalence  of  the  two  systems  (12)  and  (13), 
the  adaptive  law  for  (12)  Is  generated  In  an  identical  manner  to  that  of  (13) 
and  la  given  by  equation  (8) 

e - -re1(t)v  - -re1(t)v  r - rT  > o (14) 

Hence,  a minimal  order  adaptive  observer  can  be  designed  by  generating  suitable 
vectors  v(t)  and  w(t),  where  v(t)  is  used  in  the  adaptive  law  and  the  vector  w(t) 

Is  used  as  an  Input  to  the  observer.  This  in  turn  yields  the  well  known  form  of 
the  adaptive  observer  shown  in  Figure  [2]. 

Equations  (7)  and  (8)  represent  the  form  of  the  error  equations  for  a non- 
minimal  observer  while  equations  (12)  and  (14)  are  the  error  equations  for  a 
minimal  observer.  In  both  cases  the  time  derivative  of  the  parameter  error 
vector  8 is  adjusted  so  that  in  the  limit  6(t)  tends  to  zero.  In  equation  (6) 
which  la  the  error  equation  of  the  new  adaptive  observer,  no  provision  is  made 
to  directly  adjust  the  parameter  error  vectors  8^^  and  Instead,  new  adjustable 

A A 

parameter  vectors  8^(t)  and  8^(t)  are  introduced  to  estimate  8^  and  0^  and  compen- 
sate for  the  effect  on  e^(t) . It  is  this  fact  that  distinguishes  the  new  adaptive 
observer  from  those  that  are  currently  In  existence. 

The  error  equations  which  describe  the  behavior  of  the  new  observer  are  of 
the  form: 

e - Aae  + e1xlp  + e2u  - d[ejv  + ejv*] 

e 

6l  " (15) 

e2  " f2IV*lp’U’vl'v21 


Figure  2 The  Adaptive  Observer 


Nggm 
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where  f^  and  have  to  be  suitably  chosen  so  that  e^(t)  0,  0j(t)  -*■  and 

A 

e2(t)  0 as  t ♦ •*. 

4.  Mathematical  Preliminaries: 

The  following  lemmas  are  found  to  be  directly  related  to  the  questions 
raised  at  the  end  of  the  last  section  and  hence  to  the  final  structure  of  the 
adaptive  observer. 

1:  A time- Invariant  dynamical  system 


z ■ Az  + 6u 


Zl(t)  - h z(t) 


(16) 


where  A is  an  (nxn)  stable  matrix  In  output  canonical  form,  h - [1,0,. ..,0], 

♦T  • ♦ ] is  zero  state  Input-output  equivalent  to  the  dynamical  system 

1 2 n , 

of  order  (2n-l) 


• T 1 
w ■ Aw  + d$  v 

• 

v • Dv  + bu(t) 

1 T ' U 

v - L 


w^(t)  ■ h w 


(17) 


where  the  (n-1  x n-1)  matrix  D is  stable,  v is  an  (n-1)  dimensional  vector, 
bT  ■ [0,0,. .. ,0,1]  dT  ■ [l.dj.d^ dQ]  and  the  matrices  D and  L are  defined  as 


0 10. 
0 0 10 
0 0.. 
■1  -d-  , , 


* "d- 


L - 


1 -d2  -d3 


• • 


1 

0 


-d 

r 

• 

0 

1 


Proof : The  transfer  function  of  (16)  Is 


4 
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and  the  transfer  function  of  the  system  (17)  is 


T2(s) 


n-l  , n-2 

s + d_s  + . . .+  d 
2 n 

sn  + a, sn  *>. . .+  a 
1 n 


, n-l . , , 

f.s  +...+$ 

1 n_ 

n-l.  . . 
s +. . .+  d 


Hence,  for  zero  Initial  states  the  two  systems  have  the  same  output  corresponding 
to  any  input  u(t). 

In  the  dynamical  system  described  by  the  differential  equations 


z - Ax  + $u(t)  - d[*T(t)v1(t) ) 

a 

$ - Tz1(t)v1(t)  r - rT  > 0 (18) 

z2(t)  - hTz(t) 

where  the  matrix  A and  the  vectors  $,d,  and  v^(t)  are  as  defined  in  lemma  1, 

T -1 

r is  any  symmetric  positive  definite  matrix  and  h (sI-A)  d is  strictly  positive 
real,  4(t)  $ if  the  input  is  sufficiently  rich  and  the  overall  system  is 

asymptotically  stable. 

When  the  input  u(t)  = 0 equation  (18)  has  the  same  form  as  equation  (7)  and 

A 

hence  lk(t)D  is  bounded. 

Considering  the  differential  equations 

Ci  “ A£x  + ♦u(t)  ^ - hT51  (19) 

end  C2  - A£2  - d[iT(t)v1(t))  C2  " h\  (20) 

we  have  ^ + C2  * z and  ‘’l  + ^2  " ln  equation  (18). 

Equation  (19)  is  input -output  equivalent  to 

C3  - AC3  + d^T(t)v1(t)  C3  - hTC3  (21) 

Hence  equation  (18)  has  the  ssa te  input-output  behavior  (for  zero  initial  conditions) 

as 
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n - An  + d[(*-^(t)}Tv1(t) J 
T 

Hj  ■ h n ■ *j(t). 


(22) 


Hence,  by  the  prototype  discussed  in  section  (3)  if 

$ ■ re^tjv^t)  (23) 

and  the  input  u(t)  is  sufficiently  rich 

tim  $(t)  - $(t)  and  lim  e (t)  ■ 0. 
t-*»  tr*»  1 

5.  The  Adaptive  Observer: 

The  results  of  section  4 can  be  directly  applied  to  the  adaptive  observer 
error  equations  given  in  (15) . In  such  a case  we  have 

e ■ A e + 0,x  + 6 u - d^v^+e^v2] 

m 1 lp  2 1 2 

WT 

e,  ■ h e 

1 (24> 

- r1e1(t)v1(t) 

a ~ 

02  - r2e1(t)v2(t)  Fi  “ > 0 1 * 1.2 

1 2 

It  is  obvious  that  v (t)  and  v (t)  can  be  generated  using  two  Identical  sys- 

1 2 T 

terns  with  x,  (t)  and  u(t)  as  the  Inputs.  If  the  2n  components  of  [v  (t) ,v  (t)] 

IP 

are  linearly  Independent  and  the  input  u(t)  is  sufficiently  rich,  the  parameters 

A A 

0,(t)  and  ©2 ( t)  converge  to  the  true  parameter  errors  0^^  and 

From  the  arguments  of  section  4,  it  is  clear  that  as  t -*■  ®,  e^(t)  -*-0. 
However,  the  other  components  of  the  error  vector  e(t)  do  not  tend  to  zero  indi- 
cating that  x (t)  is  not  an  estimate  of  x (t) . To  obtain  such  an  estimate  x (t) 
m p P 

of  x (t) , the  estimates  of  the  parameter  errors  are  used  in  equation  (4)  to  yield 
P 
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1 


ft  (t)  - A ft  + bu  + 0,x,  + 0 _u 

p m p m 1 lp  2 


(25) 


The  complete  description  of  the  observer  is  then  given  by  the  differential 
equations 

• “T  1 'T  9 

(26a) 


• *T  ] 'T  2 

x * A x + b u + d[0.v  + 0,v  ) 
m mm  m 1 i 


x . - h x 
ml  m 

• A a 

x »Aft  + b u + 0,x,  + 0«u 
p m p m 1 lp  2 


(26b) 


01  - -1^ 


8 2 “ "r2elV 


(26c) 


The  resulting  observer  structure  is  shown  in  figure  3. 

Comments  on  the  Mew  Adaptive  Observer; 

a)  As  seen  from  the  lemmas  in  section  k,  the  approach  used  in  the  design  of  the 
new  observer  is  to  represent  the  unknown  transfer  function  as  the  product  of 
two  transfer  functions,  one  of  which  is  stable  and  the  other  strictly  positive 
real.  Using  such  a representation t the  prototype  error  model  developed  for 
the  design  of  non-mlnlmal  observers  can  be  directly  used  to  estimate  the 
unknown  parameters. 

b)  The  estimate  of  the  state  of  the  plant  is  obtained  by  using  the  estimates 

of  the  unknown  parameters  in  equation  26b.  The  error  in  the  state  estimate 
A 

*(t)  - x (t)  - x (t)  satisfies  the  differential  equation 
P P 

e - A^e  + l»1-61(t))x1(t)  + 

(e2-e2(t)]u(t) 

and  tends  to  zero  as  the  parameter  errors  tend  to  zero. 


GENERATOR 
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c)  The  examples  discussed  In  the  next  section  Indicate  that  the  adaptive  observer 

converges  rapidly  when  the  system  Is  of  low  order  or  when  either  the  poles  or 

the  zeros  of  the  plant  transfer  function  are  known.  However,  when  all  the 

parameters  of  a high  order  plant  are  unknown  the  rate  of  convergence  Is  found 

to  be  rather  slow  which  may  be  attributed  to  the  fact  that  the  components  of 
1 2 

the  vector  v (t)  and  v (t)  are  linearly  dependent.  Faster  schemes  of  con- 
vergence are  currently  being  Investigated. 

1 2 

d)  The  n-vectors  v (t)  and  v (t)  can  be  generated  using  identical  dynamical  sys- 
tems with  x^Ct)  and  u(t)  as  Inputs.  This  Implies  that  a single  vector  d has 

T -1 

to  be  chosen  such  that  h (sI-A^)  d is  strictly  positive  real.  In  practice, 

for  faster  convergence.  It  may  be  desirable  to  use  a correction  signal  of 

the  form  dj^v*  + ^202v^*  d[6^v^  + e2v^'  In  cflse  tbe  tran8“ 

T —1  T —1 

fer  functions  h (sI-A  ) d.  and  h (al-A  ) d0  must  be  positive  real. 

ml  m l 

6.  Examples : 

Several  plants  were  identified  using  the  new  scheme  on  a digital  computer. 

Three  typical  cases  are  Included  below.  In  the  first  case,  the  plant  Is  of 
second  order  and  has  four  unknown  parameters.  In  the  second  and  third  cases, 
the  plant  is  of  third  and  fourth  order  respectively.  In  the  first  case  the  transfer 
function  of  the  plant  Is  shown  and  the  evolution  of  the  unknown  parameters  to  their 
final  values  is  indicated.  The  states  of  the  plant  and  the  corresponding  estimates 
are  also  tabulated.  For  the  other  two  cases  the  final  values  of  the  parameters  are 


given. 
Example  1: 

w(s) [Transfer  function  of  the  Plant] 

x 2(8+2) 

2 

s +3.58+2.5 

Time 

Plant  Parameters 

in  Seconds 

Sj^  " 3.5 

a2  - 2.5  bl  " 2 

b2  - 4 

OBSERVER  PARAMETERS 

2 

3.06 

2.97  2.90 

3.99 

4 

3.53 

2.28  1.86 

3.99 

6 

3.51 

2.35  2.02 

3.98 
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Th«  following  cable  shows  the  states  of  the  plant  and  the  corresponding  estimates 
at  t • 2,4  and  6 seconds. 


Time 

X, 

lp 

lp 

X, 

2p 

2p 

2 

1.69 

2.12 

3.69 

5.65 

4 

- .14 

- .13 

- .47 

- .14 

6 

- .87 

- .87 

-1.60 

-1.47 

Kr ample  2:  In  this  case  the  plant  Is  of  third  order  with  a transfer  function 


w(8)  - ^±5g±6) 

s3+8.3s2+21.4s+16.8 

The  unknown  parameters  are  and  a^,a2  an<^  a3*  The  parameter  estimates  given  by 
the  observer  after  twelve  seconds  are  shown  below: 

A 

b^  ■ 2 bj^  - 1.9  ; Sj^  - 8.3  ^ - 8.02  ; 

a2  - 21.4  &2  “ 21.6  ; a3  - 16.8  I3  - 17.2  . 

Example  3:  A fourth  order  plant  has  a transfer  function 

w(8)  . 20[s3+2s2+2s+1] 

s +12s  +45e>60s+20 

where  the  zeros  of  the  plant  are  known.  The  following  estimates  of  the  unknown 
parameters  were  obtained  after  6 seconds. 

ax  - 12  8X  - 11.98  ; a2  - 45  t2  - 45.43 
a3  - 60  I3  - 60.39  and  aA  - 20  - 20.19 

In  all  cases  a square  wave  Input  of  amplitude  5 and  frequency  5 radlana/sec  was 
used. 

Conclusion: 

A new  structure  for  an  adaptive  obeerver  le  proposed  In  this  paper.  The 
unknown  pareaeters  of  the  eystam  are  continuously  estimated  in  a stable  feedback 
loop  and  are  In  turn  used  to  generate  the  estimate  of  the  state  of  the  plant. 


I 
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Th«  overall  system  Is  found  to  be  no  more  complex  than  the  non-minimal  observer 
currently  in  use. 

When  the  plant  can  be  represented  by  a differential  equation  of  low  order 
or  when  only  the  poles  or  zeros  of  Its  transfer  function  are  unknown  the  observer 
parameters  are  found  to  converge  very  rapidly  to  their  final  values.  For  high 
order  plants  whose  poles  and  zeros  are  not  known  convergence  appears  to  depend 
critically  on  the  auxiliary  signals  used  in  the  adaptive  laws. 

The  importance  of  the  new  adaptive  observer  lies  in  its  potential  application 


I 
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